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In this manuscript we reply to the criticisms of our paper [Eur. Phys. J. A 24, 
437 (2005)] raised in a recent preprint by Achasov and Kiselev hep-ph/0606268 and 
demonstrate that all their criticisms are completely irrelevant and misleading. 
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In a recent paper [1| we considered the radiative decay <p — > yao/Zo m the molecular 
(KK) model of the scalar mesons (ao(980), /o(980)). In particular, we showed that there 

Ph. 

was no considerable suppression of the decay amplitude due to the molecular nature of the 
scalar mesons. In addition, as a more general result we demonstrated that, as soon as the 
vertex function of the scalar meson is treated properly, the corresponding decay rates become 
very similar to those for pointlike (quarkonia) scalar mesons, provided reasonable values are 
chosen for the range of the interaction. We also confirmed the range of order of 10~ 3 -j- 10~ 4 
for the br ancMn g r atio obtained b Refs. QSQ within the m0 ,ee U . ar .node,. 

Our results imply that differences in the concrete calculations of the rate for — > yao/Zb 
for a molecular and a quark (qq or qqqq) structure of the light scalar mesons can be only 
hidden in the effective coupling constants. In either case — regardless whether the scalars 
are molecules or not — the algebra how to calculate the rates is identical. Furthermore, the 
form factors representing the involved (KK) vertices should be similar. In fact, we have 
argued that the effective couplings extracted from the <fi radiative decays seem to be more 
consistent with those expected for a molecular state or, at least, for a state containing a large 
admixture of the KK molecule. This is not surprising since for the kaon loop mechanism of 
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FIG. 1: Diagrams contributing to the amplitude of the radiative decay (j) — > 7ao//o- 



radiative decay to be operative, one needs the scalar effective coupling to be sizeable, which, 
in turn, requires a sizeable admixture of the KK in the scalar wavefunction, as demonstrated 
in P|. In other words, it means that the scalar radiative decay in the kaon loop mechanism 
proceeds via the molecular component of the wavefunction. 

As a reaction to our work a critical comment appeared (], where the authors dispute 
our results and, specifically, where they claim that our paper l| is "misleading". We shall 



refute the arguments of the authors of Ref. p] below in detail one by one. 

But first, let us briefly recall the essentials of the formalism. If the scalar meson is 
considered to be point-like, the current describing the transition — > (S = ao or / ) is 
given by the diagrams (a)-(c) of Fig. 1, and can be written as 

eg^gs 



M,, 
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I(a,b)[e v {p ■ q) -p v (q ■ e)], 



(1) 



where p and q are the momenta of the meson and the photon, respectively, is the kaon 
mass, gtj, and gs are the <pK + K~ and SK + K~ coupling constants, e u is the polarisation 
four-vector of the <t> meson, a = and b = ^f- (we deal here with stable scalars — see 
the comment on nonstable scalars below). The amplitude is transverse, M v q v = 0, and 
it is proportional to the photon momentum. The form (0) is well known, and the details 
can be found, for example, in Refs. . Specifically, the total matrix element is 

finite if evaluated consistently with gauge invariance and the integral I (a, b) entering Eq. (0) 
remains finite in the limit a — > b. 

Now, if the scalar S is treated as a non-point-like object, a corresponding form factor is 
to be attached to the K + K~S vertex. Note that the appearance of non-point- like vertices is 
in full agreement with the requirements of Quantum Field Theory (QFT). It can be derived 
naturally as a solution to the Bethe-Salpeter equation in any suitable QFT-inspired model 



3 



for the interaction of kaons which is able to bind a KK molecule. Thus it does not, by 
any means, stem from ". . . wrong interpretation of the result of calculation in the point like 
case. . .", as was claimed in Ref. 0. Also an extra graph (graph (d) of Fig. 1) should be 
added, which accounts for the additional flow of the charge. The form for this extra contact 
vertex was suggested in Re, Q, and is based on the ndnima.1 substitute considerations. 
Indeed, let the scalar vertex r depend on the total momentum of the scalar P and on the 
momentum k of the K + . Then, making the replacement —>■ — eA^, one obtains the 
"interaction current" 

-ja&n. ( 2 ) 

On the other hand, the authors of the preprint jfj] repeatedly claim that gauge invariance 
of the matrix element should be imposed by hand, namely by a subtraction at the point 



q = 0. For example, in the paper 



12j, devoted to the calculation of the <fi — > jS width in 



the molecular model, it is written explicitly: ". . .the decay amplitude 

T(p,q) = M(p,q)-M(p,0), 

M(p, q) = Mt(p, q) + M 2 (p, q) + M 3 (p, q), (3) 

where p and q are the four- momenta of the meson and photon respectively. ... The 
amplitudes Mi, M 2 and M 3 correspond to the diagrams Fig. l(a)-(c), respectively. . .." 

Let us start by stressing that, generally speaking, both prescriptions (J2J) and (jHJ) fail 
to restore gauge invariance of the amplitude. Indeed, the sum of the graphs depicted at 
Fig. l(a)-(c) is given by M v = eg^gsE^J^ with 

•V = / -^(2k + 2q-p)^2k + q) 1/ r(k,p-q)S(k)S(k-p + q)S(k + q) 

0^(2k - p),(2k -2p + q)„r(k,p- q)S(k)S(k -p + q)S(k - p) (4) 

J4 



2g^J j^r(k,p-q)S(k)S(k-p + q), 



where 



S(k) = - - (5) 

k 2 -m 2 +i0 K J 

is the propagator of the K meson, where here and in the following we use the notation 
m = rriK for simplicity. Gauge invariance requires that the amplitude is transverse, while 

QuJ^ = J -0^S(k)S(k-p)(2k-p)^[r(k + q,p-q)-r(k,p-q)]^O, (6) 
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and both prescriptions (J2J) and ® are obviously unable to repair this. 

But also here the solution to this problem is well-known in the literature (see, Ref. jl^ . 
where we list a few of the papers devoted to this subject). In principle, one is to take 
into account the strong interaction of the KK pair in the intermediate state, i.e. with the 
interaction kernel responsible for the bound scalar state formation. Inserting the photon line 
inside these graphs (the so-called structure emission) restores gauge invariance. However, in 
the soft-photon limit, one can safely do without such complications. Indeed, note that for 
soft photons (which is obviously the case for the radiative decay under consideration), one 
may safely use the current (J2J), so that, in the prescription (J2J), the piece that violates gauge 
invariance in Eq. (jSJ is cancelled by an extra contribution to J^ v , 

coming from the contact vertex (J2J) (through the contact gra ph depicted in Fig. 1(d)). There- 
fore, in the prescription (dp, the total current is transverse |l4j in the soft-photon limit, i.e. 
up to corrections of the order u/m with u being the photon energy. 

Consider now the prescription (J3J). Gauge invariance indeed requires the amplitude with 
neutral particles in the initial and final state to vanish at q — 0. However, this alone 
does not at all suffice to ensure gauge invariance, as seen from (jHJ). In addition such a 
subtraction procedure is not unique. The subtracted piece is only defined up to g-dependent 
contributions that vanish in the limit q = 0, but which might spoil the transversality of the 
transition current. Thus, there are corrections in the order of u/m too. And, finally, the 
properly calculated amplitude must vanish at q = by itself, without extra subtraction, finite 
or otherwise . Thus, based on those considerations one can also put into question the entire 
subtraction procedure suggested and advocated by one of the authors of the preprint 
However, for the present discussion it is important to note that, in the soft photon limit, 
the subtraction prescription works and, as a matter of facts, it provides results identical to 
those given by the prescription To show this, let us blindly use the prescription Q. 

. Q 

lr The case of light-by-light scattering mentioned in jy] cannot be considered as a counter-example: the 
amplitude calculated in the dimensional regularisation scheme or in the background field method nicely 
vanishes automatically in the limit uj — > 0. 
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The expression (Jl| for the transition current reads at q = 0: 



J^(g = 0) = / ^(2A;-p) M (2A;) i ,r(A;,p)5(A ; )5(A ; -p)5(A;) 



bund for the other prescription. There- 
regarding the alleged violation of gauge 



+ J ^(2k-pU2k-2p) u r(k,p)S(k)S(k-p)S(k-p) (8) 
-2^y ^r(fc,p)5(A;)5(A;-p). 

If the Ward identity, 

2k v S(k)S(k) = (9) 
is used in the first term above and an integration by parts is performed then the result reads 

J„{q = 0) = -J ^ S i k ) S i k -P)i 2k ~P)^^^ ( 10 ) 

which is obviously identical to the expression (JJJ) 
fore, we conclude that the criticism made in Ref. £ 
invariance by the amplitude derived with the help of the prescription (J2J) is rather question- 
able. 

After this lengthy introduction let us discuss the main formula of our paper [lj]. It is 
argued there (and confirmed by actual calculations) that, since the amplitude is finite even 
for the point-like limit, the range of convergence of the integrals involved is defined only by 
the kinematics of the problem. In particular, if both masses, i.e. that of the vector and of 
the scalar mesons, are close to the KK threshold, the integrals converge at ko ~ m and for 
nonrelativistic values of the three-dimensional loop momentum k, \k\ <C m. This allows us 
to perform a nonrelativistic reduction of the amplitude in the rest-frame of the (f) meson. 
The contact scalar vertex is taken in the form (J7j), yielding for the individual graphs of 
Fig. 1: 

(a ) _ (6 ) _ i f d 3 k ki{k-\q) k r{\k-q/2\) 



1 ' m3j W 3 [E v -% + iO)[E s -^ + iO) 

d 3 k r(k) 

' ° m 

d 3 k hk k 1 dr{k) 



4 ] - 2m2 tik] {2n) s Es _V +t0 > « n ) 



j(d) = % _ 

ik ' 2m 2 J (2vr)3 e v - — + iQk dk ' 
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where Ey = my — 2m, Eg = m$ — 2m. Integration by parts was performed for the last 
integral, yielding 

{d) i s f d 3 k r{k) i s f d 3 k k 2 r{k) 



Jik - 2m2 ?>ik J (27r)3 Ey _ g + . Q + 3m3 ^ J (27r) 3 (Ey-% + to)* ' (12) 

Surprisingly, it is this procedure of integration by parts which caused the main criticism 
in Ref. Indeed, it is difficult to understand why the authors of Ref. jgs] overlooked the 
simple fact that, by virtue of Eq. (JTUJ), even with the subtraction procedure (J3J) one would 
arrive at the very same result (|12p. Schematically, this looks as follows (for the sake of 
brevity we omit indices). Let us start with the form (JTTj) and denote the sum of the real 
parts of the integrals J*") + + J( c > in Eq. as J(q). Let us follow for a moment the 
authors of Ref. p and introduce the cut-off k to compute the integrals J(q, ho) and J^ d \k ). 
But it is important to realize that such a sharp cut-off obviously breaks gauge invariance - 
something which the authors of the preprint [6j have either overlooked or tacitly ignored. 
Then, as a consequence, the sum J(q, ko) + J^ d \ko) clearly exhibits a drastic dependence on 
the cut-off: the limiting value J(q, oo) + j( d \oo) is not reproduced with the values of ko up 
to several GeV (see Figs. 2, 3 of Ref. 0). This behaviour is demonstrated in Ref. [f| and it 
is advertised as being the collapse of the approach suggested in Ref. [1]. However, one can 
easily see that the sum J(q, oo) + J^(oo) vanishes at q = 0, while the sum J(q, k ) + j( d >(ko) 
does not if is taken in the form (jlljl . Thus, as already mentioned before, such a cut-off 
imposed by brute force does not respect gauge invariance and consequently, and in line with 



the arguments of one of the authors of the preprint 



calls for the subtraction: 



J R (q, k ) = J(q, k ) + J {d \k ) - [J(q, k ) + J( d ){k )}\ q=0 = J(q, k ) - J(0, ko). (13) 

If then is chosen as defined in Eq. (JT2J) one obtains J(0, ko) = J^ d \k ) and, thus, 
the procedure of integration by parts after performing nonrelativistic reduction is perfectly 
fustified even in the subtraction scheme © advocated by one of the authors of the preprint 
jy]. The regularised integral J R (q,k ) converges rapidly and for low virtual momenta of 
the kaons - which follows both from our findings and from the one of the preprint p] (see 
Fig. 4 of Ref. Jf|). In this context, let us note that the authors of Ref. 0] have deliberately 
forgotten to mention that, for realistic values for the range of the interaction, /3, i.e. (3 ~ 0.6 
GeV or larger, the result for the matrix element is actually very close to the one of the 
point-like theory. 
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Therefore, we conclude that there are no flaws in the derivation presented in Ref. 111. 
This is in a drastic contrast to the statement of Ref. |6J , where it is claimed that high kaon 
virtualities are dominating the radiative transition amplitude in the point-like limit. Note 
that, if this would be indeed the case, then the kaon loop mechanism for the radiative decay 
must be completely irrelevant there are no infinitely compact states in nature. Besides, the 
convergence of an integral for values of the loop momentum much larger than all the natural 
scales present in the problem looks rather suspicious from a physical point of view (and is 
also aesthetically unacceptable). 

A further claim by the authors of Ref. la] is that "... there is another sloppily build 

n n 

place in Ref. The authors of Ref. |l( calculate A((f)(p — > / yS(p') at q = 0, p 2 = m 2 ^ and 
{p') 2 = -nig {ms — 2m < 0), that is, at p ^ p' + q. . . ." (see p. 7 (reference [6]) of Ref. jf|). This 
is simply not true and it remains completely unclear to us what made them to draw that 
conclusion. The amplitude was calculated retaining the full dependence on q in Eq. (fTTj). It 
turned out, however, that, due to the inequality 

m{E v - E s ) > q 2 « {E v - E s ) 2 , (14) 

one may safely omit numerically the terms containing q in the expression (fTTj) . This does 
not mean that the calculations were performed under the assumption p ^ p' + q. 2 

Yet another wrong statement of Ref. (61 concerns the problem of unstable scalars. It 
remains unclear how it is related to the "question of principle" raised at p. 7 of Ref. |6J, 
but the actual situation is as follows. An attentive and unbiased reader would have noticed 
that, while in the main body of our paper [l[ the case of the stable scalar with mg — 2m < 
is considered, the Appendix is entirely devoted to the case of the unstable scalar which can 
decay into a pair of light pseudoscalars P\Pi. Contrary to the claims of Ref. jf|, the P1P2 
invariant mass distribution is calculated, both in the differential form and also integrated 
ove, the nevoid maS8 regi o„, 900 MeV< M < ^ The r esu lt8 preyed in Q 
appear to be robust against the inclusion of the finite width of the scalar. 

To summarise, it is the statement of the preprint \^ to be completely misleading. There 
are no reasons at all to believe that the radiative decays of the <j) meson into scalars would 



2 Strictly speaking, it is not legitimate to keep subleading terms in q, as both prescriptions J5J) and 10 
maintain gauge invariance only in the soft photon limit, as was discussed above. 
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allow to discriminate between molecules and compact states. 

Finally, let us reiterate the main results Ref. jl|. There, we have shown that there is no 
appreciable suppression of the <p jS width in the molecular model for the scalar mesons, 
provided reasonable values are chosen for the range of the interaction, f3 ~ 0.6 — 0.8 GeV. 
In the molecular model this range is defined by the mass of the lightest meson participating 
in the exchange force, that is by the p meson. In recent kaon-loop calculations (3, 0, Q] 
where the scalars are considered as dynamically generated states, that is, as molecules, a 
range of the interaction of similar order was used, and a good description of the data on 
the radiative transitions was achieved. If the scalars are treated as qq states, the range of 
interaction can be estimated from the range of the vertex responsible for the strong decays 
of the mesons. This can be deduced, for example, within the 3 Po model for pair creation [15j 
to be f3 ~ 0.8 GeV. Therefore, it is not possible to distinguish between a molecular or a qq 
nature of the ao(980) and /o(980) from data on radiative decays based on the value of the 
relevant interaction range because they are essentially the same. The only real difference 
in the concrete calculations of the rate for <p —>■ jcto/fa for a molecular and a quark (qq 
or qqqq) structure of the light scalar mesons is hidden in the effective coupling constants. 
The "compactness" advocated in Ref. |6j as a criterion for the discrimination is completely 
irrelevant for the radiative decay. 
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